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ABSTRACT 

The lepton-sepcific (or type X) 2HDM (L2HDM) is an attractive new physics candidate 
explaining the muon g — 2 anomaly requiring a light CP-odd boson A and large tan (I. This 
scenario leads to r-rich signatures, such as 3r, 4r and 4r -|- W/Z, which can be readily 
accessible at the LHC. We hrst study the whole L2HDM parameter space to identify allowed 
regions of extra Higgs boson masses as well as two couplings Xhaa and which determine 
the 125 GeV Higgs boson decays h —)■ and h -A- AA/AA*{t^t~), respectively. This 

motivates us to set up two regions of interest: (A) mn ^ and (B) niA ~ ~ 

C(100)GeV rnn, for which derive the current constraints by adopting the chargino- 
neutralino search at the LHC8, and then analyze the LHC14 prospects by implementing 
r-tagging algorithm. A correlated study of the upcoming precision determination of the 125 
GeV Higgs boson decay properties as well as the observation of multi-tan events at the next 
runs of LHC will be able to shed light on the L2HDM option for the muon g — 2. 
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I. INTRODUCTION 


The muon g — 2 anomaly has been a long standing puzzle since the announcement by 
the E821 experiment in 2001 [1]. During the past 15 years, development in both experimental 
and theoretical sides has been made to reduce the uncertainties by a factor of two or so, 
and thus establish a consistent 3a discrepancy. Although not signihcant enough, it could 
be a sign of new physics beyond the Standard Model (SM). Since the hrst announcement 
of the muon g — 2 anomaly, quite a few studies have been made in the context of two- 
Higgs-doublets models (2HDMs) [2-5]. Recently, it was realized that the “lepton-specihc” 
(or “type X”) 2HDM (L2HDM)* with a light CP-odd Higgs boson A and large tan/? is 
a promising candidate accommodating a large muon g — 2 while escaping all the existing 
theoretical and experimental constraints [7]. Some of the following studies showed that the 
allowed parameter space is further resrticted, in particular, by the consideration of the 125 
GeV Higgs boson decay to light CP-odd Higgs bosons h —)■ AA [8], and the tau decay 
r —)■ fiuu combined with the lepton universality conditions [9]. 

In this paper, we attempt to make a thorough study of the whole L2HDM parameter 
space in favor of the muon g — 2 explanation, and analyze the LHC tests of the favoured pa¬ 
rameter space leading to r-rich signatures like 3r, 4r and Ar + WjZ. First, we show how the 
SM Higgs exotic decays h —)■ AA as well as h —)■ AA*{t^t~) constrain the parameter space. 
It is connected to the determination of the allowed ranges of the normalized tau (lepton) 
Yukawa coupling in the right- or wrong-sign domain, and thus more precise measurement 
of the 125 GeV Higgs boson properties will put stronger bounds on the L2HDM parameter 
space. As we will see, the hAA coupling can be made arbitrarily small by a cancellation for 
mu 3> ruA only in the wrong-sign limit of the tau Yukawa coupling [10], and it opens up the 
region of niA < rnh/2 [8]. In the region of niA > ^71^/2, the three-body decay h —)■ At^t~ 
should be suppressed and the SM (right-sign) limit of the tau Yukawa coupling is allowed for 
rriA ^ 70 GeV. The allowed parameter space is further restircited by the lepton universality 
tests of HFAG which measures the leptonic decay processes at the level of 0.1% [11]. For 
this, we improve the analysis of [9] to single out proper constraints on the tree and loop 
contributions to the tau decay. 

After scanning the L2HDM parameter space, we identify two allowed regions: (A) 
the well-known region of uia mn ~ mH± and (B) uia ~ ~ (!l(100)GeV uih- 

Most of these parameter regions predict r-rich signatures easily accessible at the LHG, and 
thus can be readily probed. As a hrst step, we investigate how the current LHG 8 TeV 
data constrain the two regions, and show that the most stringent constraint comes from the 
chargino-neutralino searches. We found that the region (B) has already been excluded at 
95 % GL. For the region (A), most of the allowed L2HDM parameter region can be probed 
soon at the next runs of LHG. 

The paper is organized as follows. In Section H, we introduce the L2HDM to provide 


* In the scale invariant 2HDM with one Higgs doublet triggering electroweak symmetry breaking, the heavy 
Higgs bosons should be around 400 GeV [6], which is excluded in the type-H but not in the type-X. 
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useful formulas, and explain why a large {g — 2)^ is easily accommodated with a light CP-odd 
Higgs boson A and large tan 13. In Section III, we summarize all the relevant theoretical 
and experimental constraints, and quote some of the latest results which are not included 
in our analysis. By using the prohle likelihood method, we identify the allowed L2HDM 
parameter regions under these constraints and show them at 68% and 95% conhdence level. 
In Section IV, we discuss r-rich signatures at the LHC expected in the identihed parameter 
regions. We analyze the 3r events to identify the parameter regions excluded already by the 
current LHC 8 TeV data. In addition, we show the prospect for the future LHC14 run with 
a dedicated simulation. We conclude in Section V. 


II. 2HDM WITH A LEPTON-SPECIFIC DOUBLET (L2HDM) 


Let us hrst introduce the L2HDM to present useful formulas for our analysis heavily 
relying on the paper by Gunion and Haber [12]. Among various types of 2HDMs classihed 
by the Yukawa coupling patterns of the two Higgs doublets $ 1^2 with the same SM quantum 
numbers, the L2HDM allows the following Yukawa couplings: 

-Cy = + Y%<^ieR + c.c., (1) 


where family indices have been omitted and <l >2 = io' 2 ^ 2 - This pattern may be a result of 
a discrete Z 2 [13]: $2 —)■ $2 and $1 —)■ —$1 combined with cr —)■ —cr while the other 
fermions are invariant under the Z 2 transformation. The most general form of the 2HDM 
scalar potential is given by 


IAhdm = + 77122<h2*^2 “ ml2^\^2 + h.c. -h ^Ai 

+A3 ($^$2) + A4 (^ 1 ^ 2 ) + {^As ($I<I>2)' + 




Ae (4>I$i) 


-I-A 7 ^$ 2 *^ 2 ^ ^*^[[$ 2 ^ -|-h.C.j'. 


( 2 ) 


The Z 2 symmetry enforces Ag = A 7 = 0. However, the 771^2 tenii that softly breaks Z 2 should 
be allowed. All couplings are assumed to be real. In the desired vacuum, both doublets 
acquire VEVs, denoted as Vi and V 2 for $1 and $ 2 , respectively. Large VEV hierarchy, i.e., 
tan/3 = ^ 2 /'^^! 1, is of our interest for the explanation of the muon g — 2. 

By decomposing the doublets as (uj -|- hi + iAi)/\/2)'^, we see the model has 

three mass squared matrices of Aj, Hf and hi, which can be diagonalized by two angles a 
and (3. The physical Higgs particles in mass eigenstates are given by 


A — —Sj^Ai -f- cpA2, — —spH^ + , 

h = —Sahi -p CQ,/i2, H = Cahi S(Jl2, (3) 


where Sa and sp are abbreviations for sin a and sin/3, etc. In this paper, we adopt the 
convention 0 < /3 < 7r/2 and —7r/2 < 13 — a < 'n 12. Then the SM-like Higgs boson is 
h ~ CQ,/i 2 with either positive or negative sign for c^. In the very large tan/3 limit, two 
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Higgs doublets are almost decoupled. But some degree of non-decoupling effects, encoded 
in 0 < CjS-a 1 , will play very important roles in our study. 

The mass spectrum can be calculated analytically in terms of the coupling constants 
in the Higgs potential, but practically it is more convenient to take masses as inputs and 
inversely express coupling constants with them: 


Ai 

A 2 

A 3 

A 4 

As 


rriffcl -h tan/3 

~ ^12 cot I 3 

’ 

{mjj - mDcgSg + 2mjj±SfiC^ - mj^ 
{m\ - 2m]j±)sfjCii + 

771^2 - m\spCi3 
v'^SisCf^ 


( 4 ) 


One can see that we require an intolerably large Ai ~ tan^/3m|^/n^ > (9(10^) in the large 
tan/3 region if m ^2 = 0. Thus, the soft Z 2 breaking term mi 2 needs to be non-vanishing, 
and it is determined to be m ^2 ~ '^‘h/ tan/3. The mass splittings among the extra Higgs 
bosons are controlled by two parameters A 4 ^ 5 : 

rriH^m^ + X^v^, m%+^ m\ + - \ 4 )v‘^. (5) 


Immediately, we need A 5 ~ — A 4 ~ C^(l) to get the favored mass pattern uia mn — rnH± 
by Electroweak precision test constraints. In addition, from Eq. (4) we know that in the 
large tan/3 limit we determine A 2 ~ ~ 0.26, just as in SM. 

In general, the Yukawa couplings of the hve physical Higgs bosons, h, H, A and in 
the 2 HDM are given by 

{^hJhf + inlHf - 3 ^ 075 ^ 1 /) 

H^d + + H.C. 

V V 




where / runs over all of the quarks and charged leptons, and furthermore n, d, and I refer 
to the up-type quarks (n, c, t), down-type quarks (d, s, b), and charged leptons (e, fx, r), 
respectively. Specihed to the L2HDM, we have 


CU _cd _ cl _ 

sin/3’ cos/3’ 

Ca = -^A = cot/3, = tan/3. 


( 6 ) 
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In any type of the 2HDM, the Higgs-to-gauge boson couplings read: 


ghvv = sin(/3 - a)g^y, gnvv = cos(/3 - a)g^y, gAvv = 0, (7) 

where V refers to Z and gauge bosons. For very large value of tan/3, we have — 

~ — I'CaI = tan 13, in short, the quark Yukawa couplings of H and A are 

highly suppressed while the lepton Yukawa couplings of H and A are highly enhanced. This 
feature helps to shed a light on the muon g — 2 problem while evading various experimental 
constraints. 


III. CONSTRAINTS ON L2HDM PARAMETERS 

In this section we hrst describe all the relevant theoretical and experimental constraints 
on the L2HDM parameter space. Based on these constraints we present our results in 2- 
dimensional prohle likelihood maps. The 68% (95%) contours will be presented in dark 
(light) green in all the likelihood maps. 


A. Enhanced {g — 2)^ with large tan /3 and light A 

Recent progress in determining the muon anomalous magnetic moment = {g —2)^/2 
establishes a 3a discrepancy: 

Aa^ = - af^ = +262(85) x 10“^\ (8) 

which is in a good agreement with various group’s determinations [7]. Such an excess can 
obviously be attributed to a new physics contribution. In the framework of 2HDMs, the 
Barr-Zee 2-loop correction with a light A and r running in the loop [14] can generate a large 
positive Aa^ due to an enhancement factor of in the large tan(3 limit. Let 

us note that the Barr-Zee diagram with H running in the loop gives a negative contribution 
to Aa^ and thus a heavier H is preferred to enhance Aa^. For more details, we refer the 
readers to Ref. [7]. 

We compute {g — 2)^ by using package 2HDMC [15].1 


B. Theoretical constraints 

There are several theoretical constraints; the perturbativity, vacuum stability and uni- 
tarity bounds to be considered. All of them are implemented at the weak scale. In particular, 
the hrst imposes the highest mass scale for the Higgs states. 


t Alternative option is the public Mathematica code [16]. 
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• For the perturbativity, we put the bound: |Ai| < 47 r for 

An immediate consequence of this bound can be obtained from Eq. (5): 

+ m\, (9) 

saturated for A 5 ~ —A 4 = 47 r. Assuming a small contribution from it gives the 
upper bound mH+ ~ mn ^ 900 GeV. Note that with the large tan [5 approximation, 
Ai becomes an independent parameter and its magnitude is in principle allowed to run 
within 47 r by perturbativity. 


• Vacuum stability demands 

-^1,2 >0, A3 + \/ A1A2 > 0, IA5I < A3 + A4 + \/ AiA2. (10) 


The last condition can be rewritten as A 3 + A 4 — A 5 > — ^/AiV for mu > One 
of the key features in our discussion is that the couplings and thus the upper limits 
on the heavy Higgs masses show quite different behaviors in the right-sign (SM) and 
wrong-sign limit of the normalized Yukawa coupling [10]. Using a trigonometric 
identity, can be expressed by 

~ ^jS—a tpCjS—a- (H) 

C/3 

As found at the LHC, the 125 GeV Higgs boson h is very much SM-like requiring, in 
particular, |s/ 3 _a| — 1 and |.^^| ~ 1. Notice that this can be reached in the SM limit 
tpcp-a ~ 0 (leading to the right-sign lepton coupling ^ -|- 1 ), or in the large tan/? 
limit with tpc^-a ~ 2 (leading to the wrong-sign couplig pa — 1 ). Using the relation 
( 11 ), one hnds 


, , , , _ - {sl_^ + 1 

A3 -h A4 — A5 — - - -h 


( 12 ) 


in the large tan/? limit. Now, in the right-sign limit —)■ +1), we have 


m 


H 


, -— ‘lm\ -I- m\ 

< VO .26 X 47r + —^ 


(13) 


which puts a bound mn < 250 GeV for = 0, which is consistent with [7]. On 
the other hand, in the wrong-sign limit —)■ — 1 ), mn can be arbitrarily large 

allowing a hne-tunnig -|- ~ 0. These properties will be clearly shown in 

our Figs. 2 and 3. 


• Tree-level unitarity for the scattering of Higgs bosons and the longitudinal parts of 
the EW gauge bosons. 

The numerical evaluation of the necessary and sufficient conditions for the tree- 
level unitarity in the general 2HDM has been encoded by the open-source program 
2HDMC [15]. We deal with this constraints relying on it. Here, we point out that this 
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constraint is rather loose in the following reason. In the limit of large tan P, the param¬ 
eter Ai decouples from the other parameters A 2 , 3 , 4 , 5 , and is allowed to run between 0 
and An independently. Therefore, one can always track down a value of Ai to meet the 
requirement of the tree-level unitarity without affecting any other physical observables 
signihcantly. 


C. Electroweak precision test 

Electroweak precision test (EWPT), commonly referred to as the p parameter bound, 
is taken into account by calculating the oblique parameters. S', T and U in the 2HDMC code. 
As we are interested in a splitting spectrum of A and H, H^, the custodial symmetry is 
potentially violated signihcantly. However, as analyzed in detail in Ref. [7], taking the SM 
limit s/ 3 -a —t 1, the custodial symmetry can be restored if mH± ~ for arbitrary 

value of mAipyiu) [17]. In our scan study, we reproduce the previous results as clearly 
demonstrated in Fig. 2. Let us remark that we have updated the central values, error bars, 
and correction matrix adopted in Ref. [7], using the new PDG data [18]. 


D. Light A and Higgs exotic decay 


As we are interested in the case of a light CP-odd scalar A, the SM Higgs boson can 
have an exotic decay of (i) h —)■ AA for uia < mh/2, or (ii) h —)■ AA*{t^t~) for tua > mh/2 h 
At the moment, the current LHC data on the SM Higgs boson put a strong constraint on 
the hAA coupling XhAA and On the other hand, it will be an interesting channel to test 
the hypothesis of the L2HDM explaining the muon g — 2 at the next runs of the LHC. The 
partial decay widths of these processes are 


(i) r(A ^ 44) = - iml/ml 


(ii) P(h —)■ Att) 


32n ruh 
1 A 


hAA 


mz 


1287r^ ruhv"^ 


tan^ (3 G{m\/ml), 


(14) 

(15) 


where G{x) = (x — 1) ( 2 — - log a;) -|- 


l-5a; / 

^/Ax — 1 V 


arctan 


2a; 




arctan 


X 


y/A. 


X 


The function G{x) is a very fast monotonically decreasing function with respect to x. For 
instance, we have G(0.3) ~ 0.28 to be compared with G(0.5) ~ 0.0048. 

Cenerically, Xhaa is expected to be around the weak scale hence leading to a large decay 
width at the CeV scale, which is readily excluded. To avoid this situation, one may require 


^ In type-I and type-II 2HDM, Ref. [19] studied the possibility of two-body decay mode h -A AA while the 
three-body decay mode was ignored. 
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niA > rnh/‘2^ or arrange a mild cancelation to get sufficiently small Xhaa- Interestingly, one 
can find 

XhAA ~ —{X 3 + A 4 — As)^, (16) 

where A 3 + A 4 — A 5 is given in Eq. (12). This relation says that there could be a cancellation 
among three contributions from mA,Tnh and uih- In particular, for uih 3> mh,A of our 
interest, the cancellation is obtained only in the wrong-sign limit with < —1. This can 
be explicitly seen by taking Xhaa as a free parameter (traded with Ai) and expressing the 
normalized tan (lepton) coupling as 



- vXhAA/Si3-a 


(17) 


In the limit of mn ^ rnA and Xhaa —^ 0, it can be further approximated as —m\j/{rn\j — 
m\) < —1, and thus we have < —1. ^ We demonstrate this behavior in the right panel of 

Fig. 3. 

The presence of a light A may leave hints at Higgs exotic decay through the channel 
h -A- AA{A*) —)-4r. The upper bound of the exotic branching ratio of the Higgs decay is 
known to be 60 %, however, a mildly more stringent bound on the h -A AA mode using 
multilepton searches by CMS [20] can be set: Br(/i -A AA -A At) < 20% almost independent 
on niA [21]. In this paper we impose a conservative cut Br(/i —)■ AA{A*)) < 40%. 


E. Collider and other constraints 


• Collider searches on the SM and exotic Higgs bosons 

For various Higgs constraints from LEP, Tevatron and LHC, we rely on the package 
HiggsBounds-4.2.0 [22] incorporating the most updated data on BR{h -A rr). We 
notice that the DELPHI search [23] on the process 

e+e" ^ Z* ^ AH ^ 4r, (18) 


is sensitive to our model. The Fig. 15 in the Ref. [23] shows the region vriA + ^ 

185 GeV is excluded at 95% conhdence level. 

Specihc to our study, the 125 GeV Higgs decay h -A is of particular concern as 

it can deviate signihcantly from ±1 as indicated in Eq. (17). We use the new data 
from GMS [24] and ATLAS [25], weighted by their statistic error bars: 


firr 


1.43 ±0.40 ATLAS 
0.91 ±0.28 GMS 


(19) 


^ The case with S/s-a ~ —1 (or equivalently cos a « —1), i.e., the reversed couplings of other SM couplings, 
is completely excluded from our numerical results. So, we have s,g_a ~ +1 in this paper. 



• Bs -)■ /i+yU 

The light A contribution to the decay Bg —)■ becomes sizable if uia ^ 10 GeV. 

In our analysis, we do not include this constraint as it is irrelevant for niA > 15 GeV. 
More details can be found in Refs. [8, 9]. 

• r decays and lepton universality 

In the limit of light and large tan ( 3 , the charged Higgs boson can generate signif¬ 
icant corrections to r decays at tree and 1-loop level [26]. Recent study [9] attempted 
to put a stringent bound on the charged Higgs contributions from the lepton univer¬ 
sality bounds obtained by the HFAG collaboration [11]. Given the precision at the 
level of 0.1 %, the HFAG data turned out to provide most stringent bound on the 
L2HDM parameter space in favor of the muon g — 2. Thus, it needs to be considered 
more seriously. For this, we improve the previous analysis treating the HFAG data in 
a proper way. 

From the measurements of the pure leptonic processes, r —)■ nuu, r —)■ evv and 
fx —)■ ez/z/, HFAG obtained the constraints on the three coupling ratios, {gr/gn) = 
\/T{t —)■ ez/z/)/F(yU —)■ ez/z/), etc. Dehning 5iii = {gi/gv) — 1, let us rewrite the data: 

= 0.0011 ±0.0015, = 0.0029 ±0.0015, = 0.0018 ± 0.0014 (20) 

In addition, combing the semi-hadronic processes ti/K —)■ /iz/, HFAG also provided 
the averaged constraint on {gr/g^i) which is translated into 

= 0.0001 ± 0.0014. (21) 

We will impose the above lepton universality constraints in our parmeter space. 

Now, it is important to notice that only two ratios out of three leptonic measurements 
are independent and thus they are strongly correlated as represented by the correlation 
coefficients [11]. Therefore, one combination of the three data has to be projected out. 
One can easily check that the direction — (5(.g ± has the zero best-fit value and 
the zero eigenvalue of the covariance matrix, and thus corresponds to the unphysical 
direction. Furthermore, two orthogonal directions ± (5(-g and ± 5(-g ± 25^^^ are 
found to be uncorrelated in a good approximation. In the L2HDM, the deviations 
from the SM prediction <5;;/ are calculated to be 




■'Til, '^loopi ^re ^tree “ 1 “ X^loop: 

Here 6tree and 6loop are given by [26]: 

mlmj, , ^ ^g{rri^^/ml) 


± 


Sle = S. 


tree-! 






loop- 


^tree. — 


8m 


_± R - 


^loop — ^ ^ 


8^2- 


TT 


2 


IT - {H (Xa) + (xn) + {Xk}) 


( 22 ) 


( 23 ) 
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where f{x) = 1 — 8x + 8a:^ — x"^ — 12x^ ln(a;), g{x) = 1 + 9x — 9x^ — + 6x(l + a;) ln(x), 

H{x) = ln(x)(l + x)/(l — x), and x^ = From Eqs. (20,21) and (22), one 

obtains the following three independent bounds: 


1 

71 


^tree 


+ 


loop 



0.0028 ±0.0019, 

0.0022 ±0.0017, 
0.0001 ±0.0014. 


(24) 



mjj± (GeV) 


FIG. 1: The contours of lepton universality likelihood profiled on tan/1) plane. The red, 

blue, and black lines are corresponding to 99%, 95%, and 90% confidence limit, respectively. 

Based on the constraints Eq. (24) on the two fundamental free parameters 6tree and 
Sloop, we can draw the the lepton universality likelihood contours, where we found the 
minimum value Xmin ~ 0.229. In Fig. 1, we present prohle likelihood contours on 
the mfp±-ta.ia/3 plane the red, blue, and black lines are corresponding to 99%, 95%, 
and 90% conhdence level, respectively. Note that the Sioop is always negative in the 
region of our interest listed in Table I. On the other hand, Stree depends only on the 
parameter tan /3/mH± and negative in most of the region but can be also positive. In 
a hne-tuned region located tan/3/m//± ~ 1 GeV“^ as we can see in the large tan/3 
and small mH± corner in Fig. 1, where the positive Stree and the negative Sioop cancel. 

We also found that lepton universality likelihood is practically not sensitive to the 
heavy neutral Higgs mass mu and cos(/3 — a) in our region of interest. Hence, we 
show the lepton universality contours only on the mi^±-tan/3 plane (Fig. 1) and on 
the m/i-tan/3 plane (Fig. 4 left panel). 
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Let us finally remark that we use Gaussian distribution or hard cut for the likelihood 
functions to impose the experimental constraints. When the central values, experimental 
errors and/or theoretical errors are available, Gaussian likelihood is used. Otherwise the 
hard cut is adopted such as the Higgs limits implemented in HiggsBounds. 


2HDM parameter 

Range 

Scalar Higgs mass (GeV) 
Pseudoscalar Higgs mass ( GeV) 
Charged Higgs mass ( GeV) 

C/3—a 

tan/3 

Ai 

125 < mn < 400 

10 < niA < 400 
94 < mfj± < 400 

0.0 < C/3_q, < 0.1 

10 < tan P < 150 
0.0 < Ai < dvr 


TABLE I: The scan ranges of the input parameters over which we perform the scan of L2HDM. 
Note that we adopt the convention in 2HDMC; —7r/2 < a — /3 < 7r/2 and 0 < /3 < 7r/2, and use the 
parameter Ai as an input parameter instead of 771^2 in order to make the scan more efficient. 


F. Results 

Our input parameters and the scan ranges of them are summarized in Table 1. Some 
comments are in order, (i) We focus on the case that the SM-like Higgs boson h is the lighter 
GP-even Higgs boson with mass 125 GeV [27]. ^ (ii) We require cos(q; — P) < 0.1, which 
guarantees that h couples to quarks and vector bosons without appreciable deviation from 
the SM predictions. The updated LHG results can be found in Ref. [28]. (hi) The upper 
bound on mH,H± < 400 GeV is put by hand since we are interested in the relatively light 
region testable at the LHG near future. In principle, they can be as heavy as about 900 
GeV according to the perturbativity constraints, (iv) We restrict ourselves to tan P < 150. 

We show the scan results in several 2 dimensional prohle likelihood maps from Fig. 2 
to Fig. 4. The inner green (outer light green) contours are 68% (95%) conhdence region. 
The points are summarized in the following: 

• The left panel of Fig. 2 shows two separated allowed regions. The majority is crowding 
around the line mn = rriH+i which is in well accordance with the EWPT via accidental 
degeneracy between H and H^. Note that there is a lower bond on mn ~ ^h+i about 
130 GeV. The minority is on the small island with quite light near mH± ~ 100 
GeV, just in the vicinity of the LEP bound on charged particles. With the help 


^ We have checked the case that the SM-like Higgs is the heavier CP-even Higgs. We found that the allowed 
region is rather restricted at mh — mn — 125 GeV, which is the similar solution to the subset of region 
(B). 
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of the right panel of Fig. 2, one hnds a mild degeneracy between A and II^ with 
niA ~ 100 — 180 GeV and mff± < 160 GeV. For tua > 100 GeV, tan (3 needs to be 
larger than about 70, see Fig. 4. We call the former region as Region A and the latter 
as Region B. Note that the fragmentation of the plots, particnlarly in the region B of 
the left panel of Fig. 2, is dne to a coarse-tnning likelihood. As we will see in the next 
section. Region B is already exclnded by the cnrrent LHG 8TeV data. 

• The left panel of Fig. 3 shows the relation between XhAA and We see only |A/iAyi| ~ 

0 is allowed for tua < 60 GeV, while larger is allowed for rriA ^ 60 GeV. The 

right panel of Fig. 3 shows the relation between vs. rriA- In the region ttia ^ 70 GeV, 
only the wrong-sign region < 0) is allowed. It is consistent with snppressed XhAA 
seen in the left panel as discnssed in Eq. (17). For heavier A, there appears the 
right-sign region. 

• Remarkably, the < 60 GeV region tends to show an enhancement in Br(h —)• rr), 
np to a factor ~ 4. While above it both (mild) enhancement and snppression 
are possible. Farther precise measnrement of Br(h —)■ rr) helps to shrink the allowed 
parameter regions. 

• In the left panel of Fig. 4, The contonrs of lepton nniversality likelihood are also 
presented in 99% (red), 95% (bine), and 90% (black) conhdence limit. The region 
with tan f3 < 140 with small tua allowed by other constraints are very constrained 
by lepton nniversality. However, the region located at the large tan (3 > 140 are 
always allowed by the hne-tnning cancellation between 5tree and 5ioop by selecting an 




FIG. 2: Features of the Higgs spectrum with a light A facing EWPT. The inner green (outer light 
green) contours are 68% (95%) conhdence region. Distribution on the mn — mf{± plane (left) and 
the niA — fn}j± plane (right). 
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(GeV) rtij^ (GeV) 

FIG. 3: The 2-dimensional profile likelihood. The inner green (outer light green) contours are 
68 % (95%) confidence region. Left panel: the coupling ^hAA (in unit of v) versus niA- Right panel: 
the reduced coupling of leptons versus ruA- 


appropriate mH±- The lower tan/? region allowed at 95% appears to be a consistent 
combination of the same 95% contour lines with different values of mH± in [9]. 

• A light A with rriA ~ 20 — 63 GeV is of our particular interest.** In this region, the 
wrong-sign limit ~ —1) has to be realized and thus the lower bound on tan/3 is 
correlated with the upper bound on cos(a — /3), which can be seen from the right panel 
of Fig. 4. We can also see that the two discrete regions correspond to the right-sign 
limit (tan/?cos(/? — a) ~ 0) and wrong-sign limit (tan/?cos(/? — a) ~ 2) as described 
around Eq. (11). 

• The exotic Higgs decay h -A- AA or h —)■ Att is a promising channel to probe the 
L2HDM explanation of the muon g — 2 as its branching ratio can be quite sizable 
unless there is a particular reason to suppress Xhaa as shown in Fig. 5. 


IV. r-RICH SIGNATURE AT LHC 


In the previous section, we identihed two favored regions of the L2HDM parameter 
space. In this section, we discuss how the current LHC search results can constrain this 
model further. Since the relationship between ttia and tan/? is constrained by the {g — 2 )^, 


** Remark again this region is further reduced by considering the tau decay and lepton universality data [9]. 
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FIG. 4: Left: Distribution on the ruA — tan/? plane (left), and the niA — cos(q; — /?) plane (right). 
The contours of lepton universality likelihood are also presented in 99% (red), 95% (blue), and 
90% (black) confidence limit. 


> 

QJ 


u 



FIG. 5: Plots of the SM-like Higgs exotic decay Br(/i —)■ AA) (for uia < fnhl2) and Br(h —>• At^t ) 
(for mhl‘2 < niA ^ ui/j.). All the scatter points satisfy the constraints described in the text in 2(t. 


as shown in the left panel of Fig. 4, we can simply parametrize tan/? as a function of 

tan/? = 1.25 (-^)+25, (25) 

which will be assumed in this section. We left with three Higgs mass parameters mAi^H^ 
and mH±, which determine phenomenologies at the LHC. 
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The bulk parameter space with tua rnu ~ is a clear prediction of the lepton- 
specihc 2HDM considered in this paper. Since the extra Higgs bosons are mainly from 
the “leptonic” Higgs doublet with a large tan/?, all the three members are expected to 
dominantly decay into the r—flavor, leading to r—rich signatures at LHC [29-31] via the 
following production and ensuing cascade decay chains: 


pp ^IF^* ^ H^A 

W 

H- 

+ 

-s 

1 

(26) 

pp -^Z*/'y* —)■ HA —)■ 

-s 

+ 

-s 

1 

W 

+ 

-s 

1 

(27) 

pp ^IF^* ^ H^H ^ 

-s 

H- 

+ 

-s 

1 

(28) 

pp ^Z*/Y H+H- 

+ 

w 

1 

(29) 


As seen in Fig. 2, we can also hnd a small island at the right-lower corner of the plot where 
~ rriA ~ 100 GeV, which we call Region B while the above bulk region we call Region 
A. In the following we £x mH± in the two regions based on the best fit point: 

Region A: mH± = ran + 15 GeV 

Region B: mH± = max(90 GeV, 0.8mA + 10 GeV) 

With these relations we explore mA-mu plane. 

A large tan /3 enhances the lepton Yukawa couplings of extra Higgses /H/A leading 
to a fast decay into tan leptons in general. The light pseudo-scalar A indeed decays into rr 
essentially at 100%, however, the heavier /H, in the presence of this light A, can sizably 
decay into AW^/Z via electroweak gauge interactions. This partial decay width is enhanced 
by the well-known factor in the limit and expressed as 

1 

for m]j+ > M^, (30) 

lOTT \ V / 

where X{l,x,y) = {1 — x — yY — Axy. It can be compared with the partial decay width of 

—)■ TV 

r(iL+ r+z/) 

lOTT \ V I 

From Eqs. (30) and (31) one can see that the WA channel turns out to dominate over the 
TV channel when mH+ > v^^^rtan/d as shown in the left panel of Fig. 6, where we plotted 
the branching ratio of —)■ AW^. We can get the decay width r{H ZA) hj replacing 

mH+ and Mw with mn and Mz, respectively, in the above expression, and its branching 
ratio is also shown in the right panel. 

Even if H/H^ undergoes the decay involving ZjW^, the associated A will eventnally 
decay into rr and thus multiple r signature up to 4r -|- IF or/and Z would be one of the 
peculiar signatures of the model at the LHG. 
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BR(H+ -» AW+) 


BR(H -» AZ) 



m/,[Gev] 


FIG. 6: Contour plot of branching ratio 
Br(fl’*’ —)• ~ 1 in Region A. The relation tan/3 



m^lGev] 

AW+) and Br(R' ^ AZ). Bi{H+ AIF+) + 
1.25771^4 + 25 is used. 




FIG. 7; 95% CL. contour from the chargino-neutlarino search at LHC 8TeV shown in uia vs. mn 
plane for Region A (left) and Region B (right). 

A. Current Constraints 

Current LHC 8 TeV data already set the constraints in the parameter space we 
are interested in. In both Region A and Region B we take model point grid with 
ruA G [20,200] GeV and mn G [140,320] GeV both with 20 GeV steps, that is, 100 model 
points for each region. We generate the 50,000 signal events with MadGraph [32] for each pa¬ 
rameter point and interfaced to CheckMATE 1.2. 0-beta [33] for checking the current bound 
with 20 fb“^ data at 8 TeV LHG. The analyses implemented in the CheckMATE are listed in 
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the Table IV A. We checked all the analyses and considered that a model point is excluded 
when at least one analysis excludes it at 95% C.L. 

Fig. 7 shows the estimated 95% C.L. exclusion contours. For most of the parameter 
space, the strongest constraint comes from the chargino-neutralino search in ATLAS [34]. 
Especially, it is from the signal region “SR2ra” therein, which requires two r leptons and 
an additional isolated lepton, with > 100 GeV, I^t > 50 GeV and &-veto. Heavier 

tuh > 200 GeV (Region A) or run > 280 GeV (Region B), and light tua < 50 GeV are still 
allowed and we will show later that the next run of LHG can explore some of the regions. 
For the heavier mn regions the sensitivities are weaker just because of the smaller cross 
sections, while for light region it is because rs from lighter A decays become softer and 
thus the acceptance quickly decreases. Moreover, the H/H^ —)■ AZjW^ decay modes also 
start open to decrease the number of hard rs from direct H/H^ decays. Note that the 
exclusion of the lighter parameter space is of interest only for Region A, since for Region 
B the interesting tua in our scenario to explain {g — 2)^ is conhned to be lie above 100 GeV 
as you can see in Fig. 2. 


B. 14 TeV prospects 

In this section we estimate the reach of the LHG 14 TeV in Region A and B based on 
the model point grids dehned previously for the LHG 8 TeV study. The signal cross sections 
depend on heavy Higgs masses, and in Fig. 8 we show the contour plots of total cross section 
on the ruA — rnu plane for Region A (Region B) in the left (center) panel. Actually, for 
relatively small tua the dominant contribution comes from the H^A production while the 
HA production contributes secondarily; HH^ and contributions are subdominant. 



FIG. 8: Total signal cross section dependence in mA vs. mu plane in Region A (left) and Region 
B (center). Right panel: Missing transverse momentum distributions for the signal benchmark 
point C {mA = 100 GeV and mn = 200 GeV in Region A) and various BG processes. 

For the Standard Model background processes we consider tt, hF+jets, Z+jets, and di¬ 
boson productions W~^W~, ZZ. All background events are generated with ALPGEN [35] 
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arXiv number 

description 

integrated luminosity [fb ^] 

atlas-1308-2631 

ATLAS, 

0 leptons -p 2 b-jets + etmiss 

20.1 

atlas-1402-7029 

ATLAS, 

3 leptons + etmiss (chargino-pneutralino) 

20.3 

atlas-1403-4853 

ATLAS, 

2 leptons -p etmiss (direct stop) 

20.3 

atlas-1403-5294 

ATLAS, 

2 leptons -p etmiss, (SUSY electroweak) 

20.3 

atlas-1403-5294-CR 

ATLAS, 

2 leptons -p etmiss CR, (SUSY electroweak) 

20.3 

atlas-1404-2500 

ATLAS, 

Same sign dilepton or 31 

20.3 

atlas-1407-0583 

ATLAS, 

1 lepton -P (b-)jets -p etmiss (stop) 

20.3 

atlas-1407-0600 

ATLAS, 

3 b-jets -p 0-1 lepton -P etmiss 

20.1 

atlas-1407-0608 

ATLAS, 

Monojet or charm jet (stop) 

20.3 

atlas-1502-01518 

ATLAS, 

Mono jet plus missing energy 

20.3 

atlas-conf-2012-104 

ATLAS, 

1 lepton -P > 4 jets -p etmiss 

5.8 

atlas-conf-2012-147 

ATLAS, 

Monojet -p etmiss 

10.0 

atlas-conf-2013-021 

ATLAS, 

WZ standard model (3 leptons -p etmiss) 

13.0 

atlas-conf-2013-024 

ATLAS, 

0 leptons -P 6 (2 b-)jets -p etmiss 

20.5 

atlas-conf-2013-031 

ATLAS; 

Higgs spin measurement (WW) 

20.7 

atlas-conf-2013-036 

ATLAS; 

4 leptons -p etmiss 

20.7 

atlas-conf-2013-047 

ATLAS, 

0 leptons -p 2-6 jets -p etmiss 

20.3 

atlas-conf-2013-049 

ATLAS, 

2 leptons -p etmiss 

20.3 

atlas-conf-2013-061 

ATLAS, 

0-1 leptons -P > 3 b-jets -p etmiss 

20.1 

atlas-conf-2013-062 

ATLAS; 

1-2 leptons -p 3-6 jets -p etmiss 

20.1 

atlas-conf-2013-089 

ATLAS, 

2 leptons (razor) 

20.3 

atlas-conf-2014-014 

ATLAS, 

2 leptons -p b-jets (stop) 

20.3 

atlas-conf-2014-033 

ATLAS, 

WW standard model measurement 

20.3 

atlas-conf-2014-056 

ATLAS, 

ttbar spin correlation measurement 

20.3 

cms-1303-2985 

CMS, alpha-T -|- b-jets 

11.7 

cms-1301-4698-WW 

CMS, WW standard model measurement 

3.5 

cms-1405-7570 

CMS, Various chargino and neutralino 

19.5 

cms-smp-12-006 

CMS, WZ standard model (3 leptons -|- etmiss) 

19.6 

cms-sus-12-019 

CMS, 2 leptons, > 2 jets -|- etmiss (dilep edge) 

19.4 

cms-sus-13-016 

CMS, OS lep 3-|- b-tags 

19.5 


TABLE II: The list of the analysis used in our analysis implemented in the CheckMATE. The list is 
found in the CheckMATE/data/ directory. 


+ Pythia [36]. We only consider leptonic decay modes including tan for all processes as 
later on we select events with at least 3 leptons including tans. To include the mis-tagging-r 
effects, we generate the MLM-matched samples [37] with 2 to 3 additional jets for Lh+jets, 
and with 1 to 2 additional jets for Z+jets. Cross sections with the above generation cut are 
102 pb, 1365 pb, 714 pb, 8.13 pb, 0.942 pb and 0.112 pb for ti, LL+jets, Z+jets, W~^W~, 
and ZZ, respectively. 

As this model predicts r-rich signatures the signal is sensitive to r-tagging and we im- 
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plement r-tagging algorithm using track and calorimeter information from Delphes 3.0 [38] 
as described in Ref. [39], which basically is a simplihed version of the ATLAS r-tagging al¬ 
gorithm [40, 41], We use two variables: 


^max = max AR(pj,Pi) and /core 

tracks 


X]r<0.2 


ZT'calo 

2l/rp 


^calo ’ 


(32) 


where pj is the jet center direction and the distance of the furthest track among pi (with 
Pt > I GeV) to Pj is denoted as i?max; is the Et deposited in each calorimeter tower 
and the summations run over the calorimeter towers within the cones centered around pj 
with cone size R < 0.1 and 0.2 for the numerator and the denominator, respectively. Both 
RmsLx and /core measure essentially how narrow the jet is; r-jet is expected to be narrow and 
gives a smaller i?max and /core ~ 1- We found these two variables are most relevant for the 
discrimination. 

We show i?max and /core distribution in Fig. 9. We also show the ROC curve obtained 
by changing the cut value for R^ax < Rmax with hxing Z™/, = 0.95 for /core > /core- 
Compared with the plot shown in Ref. [41], our simulation is reasonably conservative up 
to the signal efficiency ~ 60%. We select the working point with Rj/^x = 0.1, which gives 
er = 59% with the background jet rejection l/esc = 97. 



FIG. 9: ROC curve for our r-tagging algorithm. Our working point is denoted with a filled square, 
where 59% efficiency with 1% mis-identification efficiency for QCD jets is obtained. 

We apply the following event selection cuts to the signal and BG events. First, we 
require events with at least 3 r-tagged jets, based on the algorithm explained above. At 
this stage the dominant background becomes ti, kF-|-jets and Z-l-jets. Next, we require 
enough missing momentum > 100 GeV to efficiently reduce the kF-|-jets and Z-l-jets 
contributions. The normalized Ifr distributions are shown in the right panel of Fig. 8. 
Finally, to reduces the tt background, we veto events with any 5-tagged jet with pt > 25 GeV 
nor any jet with p^ > 50 GeV. This cut efficiently reduces the remaining backgrounds. 
Table III summarizes the number of events after the successive selection cuts in unit of fb 
for the various BG processes and for the signal benchmark model point G. We compute the 
signal to background ratio S/B and signihcance based on statistical uncertainty S/y/B. The 
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significance quoted here is based on the integrated luminosity of 25 fb“^. We can use the /i/i 
modes as suggested in Ref. [30] to improve the sensitivity and to reconstruct the events but 
we mainly focus on r-rich signatures, which require a relatively low statistics to set limit 
and expected sensitive at the early stage of LHC run 2. 

We show the results for several selected benchmark points A to F in detail. Table IV 
collects the numbers and signihcances including the other benchmark model points. 


selection cuts 

point G 

tt 

VF+jets 

Z+jets WW 

IFZ 


total BG 

S/B 

total iTgen [fb] 

153.580 

102 • 10^ 1365 • 10^ 714 • 10^ 8125 

942 

112 

2190 • 10^ 

- 

- 

ne > 3 

21.713 

273.27 

138.59 

3412.84 6.495 88.937 26.965 

3947.1 

- 

1.7 

Ur > 3 

4.386 

5.837 

13.776 

91.324 0.070 

0.343 

0.174 

111.52 

0.04 

2.1 

l^T > 100 GeV 

1.179 

1.482 

0.232 

1.244 0.000 

0.018 

0.003 

2.980 

0.4 

3.4 

Uh = rij = 0 

0.857 

0.163 

0.000 

0.505 0.000 

0.017 

0.003 

0.688 

1.2 

5.2 


TABLE III: The number of events after applying successive cut for 14 TeV LHC. Benchmark point 
C {niA = 100 GeV, ruH = 200 GeV) is shown for the signal. The significance quoted is based on 
integrated luminosity of 25 fb“^. 



point A point B point G point D point E point F 

niA [GeV] 

20 

40 

100 

40 

100 

180 

mff [GeV] 

200 

200 

200 

260 

260 

260 

total CJgen [fb] 

270.980 241.830 

153.580 

100.430 

71.271 

44.163 

n£ > 3 

6.606 

16.681 

21.713 

7.110 

11.962 

8.822 

Ur > 3 

0.894 

2.602 

4.386 

0.888 

2.346 

1.971 

> 100 GeV 

0.201 

0.547 

1.179 

0.209 

0.765 

0.926 

Uh = Tlj = 0 

0.098 

0.314 

0.857 

0.121 

0.479 

0.631 

S/B 

0.1 

0.5 

1.2 

0.2 

0.7 

0.9 


0.6 

1.9 

5.2 

0.7 

2.9 

3.8 


TABLE IV: The number of events after applying successive cut for 14 TeV LHC. The significance 
quoted is based on integrated luminosity of 25 fb“^. 

Based on the signihcance values we show the expected discovery reaches at LHC 14 
TeV in Fig. 10. The left panel corresponds to Region A and the right panel does to Region 
B. Both panels show the expected 2(T, 3a and 5a discovery reach contours with assumed 
integrated luminosity of 25 fb“^. It is seen that most of the interesting parameter regions 
can be covered. Only limitation is for the region with light rriA and heavy mn where the 
sensitivity becomes weak even though the intrinsic signal cross sections are not so small. 
The reasons are again because of the smaller acceptance for the softer r and longer decay 
chains involving Z/W as explained in the previous section on 8 TeV analysis. Moreover, in 
such a region, a light A from heavy /H decay will be boosted, resulting in a collimated 
r—pair which becomes difficult to be tagged as two separated r-jets. It is one of the reasons 
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to have less acceptance for this parameter region. We can estimate the separation Rrr of 
the T leptons from A decay: 


2m 

Rtt ^ 

Pt 


imA 


mH±/H 


1 - 2 - 


^A+^W/Z 


H±/H 


H±IH 


(33) 


For example, Rrr ~ 0.4 for mn = 300 GeV and = 30 GeV, and Rrr ~ 0.3 for mn = 
400 GeV and = 30 GeV. Since the jets are usnally dehned with R = 0.5, the r—pair 
starts overlapping. We indicated the region with the overlapping r problem in red lines in 
the left panel of Fig. 10. In that region, we have to think of how to captnre the kinematic 
featnres of the boosted A —>■ . We may be able to take the overlapping r problem 

as an advantage by ntilizing jet substructnre stndy, which is already proven usefnl [42-44]. 
For example, using di-tau tagging as proposed in Ref. [45] might be benehcial, although we 
leave this for future work. 


R„=0.5 R„=l 




FIG. 10: 2(T, 3cr and 5 cj discovery reach contours at LHC 14 TeV shown in mA vs. mn plane for 
Region A (left) and Region B (right). Assumed integrated luminosity here is 25 fb“^. Benchmark 
points selected in Table IV are indicated with filled boxes. Red lines indicate the region with 
expected smaller r separation of Rrr ~ 0.5 and 1. 


V. CONCLUSIONS 

The lepton-sepcihc (or type X) 2HDM is an interesting option for the explanation 
of the muon g — 2 anomaly which requires a light GP-odd Higg boson A and large tan/?. 
In this paper, we made a scan of the L2HDM parameter space to identify the allowed 
ranges of the extra Higgs boson masses as well as the related two couplings and Xhaa 
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of the 125 GeV Higgs boson which govern its standard and exotic decays h —)■ and 

h —)■ AA/AA*{t^t~), respectively. The tan Yukawa coupling is found to be either in the 
wrong- or right-sign limit depending on the mass of A. More precise determination of the 
standard tan Yukawa coupling and a possible observation of one of the above exotic modes 
would provide a hint for the current scenario. 

There appear two separate mass regions in favor of the muon g — 2: (A) rriA ^ 

mH± and (B) tua ~ rnH± ~ lOOGeV tuh, which lead us to set up two regions of interest 
for the LHG study: (A) mH± = m//-|-15GeV, and (B) mH± = max(90GeV, 0.8mA + 10GeV) 
with tau/d parametrized by tan(3 = 1.25(m^/GeV) -|- 25. In these parameter spaces, one 
expects to have r-rich signatures readily accessible at the LHG through the extra Higgs 
productions pp —)■ AH^/AH//HH followed by TT —)■ AZIt^t~ H^ —)■ AW^ jr^v 
and A —)■ Indeed, the current LHG8 data start to exclude (yet mild) some of the 

above two regions: rriH up to about (A) 200 GeV and (B) 280 GeV for rriA > 50 GeV from 
the consideration of the ATLAS neutralino-chargino search results. However, the region of 
i^A ^ 30 GeV (with tan f3 < 40) which also satisfies the tan decay and lepton universality 
data [9] is hardly tested by the r-rich signatures in near future even though HL-LHG should 
be able to over the region. Thus, further study, for example, on the boosted A —)■ rr will 
be required in the next runs of LHG to cover all of the L2HDM parameter space explaining 
the muon g — 2 anomaly. 
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